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We study the pseudoscalar glueball candidates in a chiral effective Lagrangian model proposed 
by 't hooft, motived by taking into account the instanton effects, which can describe not only the 
chiral symmetry breaking, but also the solution of Ua(1)- We study the parameter space allowed 
by constraints from vacuum conditions and unitary bounds. By considering two scenarios in ++ 
sector, we find that parameter space which can accommodate the ^ sector is sensitive to the 
conditions in ++ sector. From our analysis, it is found that three n states, i.e. 77(1295), 77(1405), 
77(1475), can be glueball candidates if we assume that the lightest ++ glueball has a mass 1710 
MeV. While there is no v glueball candidate found in experiments if we assume that the lightest 
++ glueball has a mass 660 MeV. 

PACS numbers: 12.39.Mk, 12.39.Fe, 11.15.Ex, 11.30.Rd 



I. INTRODUCTION 



Currently, as provided in PDG [l[ and other literature fl, H, H[, there are about 11 h particles 
identified in total, which are tabulated in Table [U It is urgent to know which are normal qq states and 
which are exotic states. For a ground qq state, in the quasiclassic picture, the quark and antiquark are 
attached to a string with the tension a. The mass squared of its higher excitation states can be described 
as m 2 = 47rcr77,. This linear n dependence can be used to classify these states properly and find out exotic 
states. Therefore, it is useful and helpful to put these states on the Regge trajectories, as shown in Fig. 
(P} where we assume that Aira — 1 which is consistent with the analysis presented in Q. We find that 
these states can be well-organized by the radial excitations of 77 and 77', except the state 77(2190), which 
is proposed as a glueball candidate in @. 

Let's comment on these resonances briefly: the first two are well-known 77 and 77'. For the following 
three, the 77(1295) is believed to be the first Regge excitation of 77 or even doesn't existjl]. The 77(1405) 
is possibly a glueball while the 77(1475) is believed to be the Regge excitations of ss [H, H|7 Or there is no 
split at all and there is just the radial excitation of 77. The 77(1760) is the radial excitations for 77, while 
the A"(1835) is assumed to be the first Regge excitation of 77' ;4j. The ?7(2070) or the 77(2190) can be 
interpreted as a pseudoscalar glueball [|| or just the second Regge excitation of 77' 0]. The last ?7(2250) 
should also be the radial excitation of 77. Interested readers can refer to @, H, 0] for further reading and 
more details. 
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TABLE I: The 11 spectra of 77 mesons (The unit is in MeV) . 



Another interesting and long-time quest is for glueball candidates. Glueball states are the consequence 
of the non-Abelian nature of QCD gauge theory. The first pseudoscalar glueball candidate was found in 
Mark II and Crystal Ball Collaboration and was dubbed as t(1440) or G(1440). Nowadays, this particle 
is named as ?7(1405) by PDG. Currently, it is widely believed that glueballs can be produced in the J/\P 
radiative decay via the "/gg channel, where the final two gluons have a considerable probability to form 
a glueball bound state. Meanwhile, a pure glueball state is forbidden to decay into a photon pair or a 
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FIG. 1: The n and r/' spectra on Regge lines. The slope of these two lines is taken as 1 GeV 2 . 



lepton pair at the leading order ( which is known as the rule of thumb OZI suppression ), therefore its 
partial decay widths to these channels should be suppressed and should be much smaller than those of a 
qq state. Such a feature can be used to distinguish hardonic and glueball states. For instance, by study 
the productions and decays of those r? states, several groups 0, Q come to the conclusion that 77(1405) is 
a good glueball candidates. At the same time, there are other states, like 77(1760) @ and A(1835) [Io| . 
which could also be a glueball state. The consensus on the question which is glueball among these h 
states has not been achieved yet. 

It is well-known a realistic glueball state is a mixture of constituent gluons and quarks. A large mixing 
prohibits us to isolate a hadronic state from a glueball state via either productions or decays. Therefore 
it seems difficult to have a definite and positive answer to the question on the r\ spectra: which spectrum 
might be the glueball candidates? Then it is urgent to have a trustable and systematic method to describe 
the glueball and meson mixing and their interaction, while the effective Lagrangian method is a good one 
for such a purpose. 

The effective Lagrangian method is a useful method for hadron physics. The advantage of this method 
lies in the fact that the global symmetries of QCD can be realized in a convenient way. It is natural to 
use this method to describe the psuedo-scalar glueball candidates. In principle, the effective Lagrangian 
should be derived from QCD, as shown by Cornwall and Soni [ll[ for pure glueball sector, or as demon- 
strated by Gasser and Leutwyler [l^, [H| for the case when only Goldstone mesons are concerned. It 
is well-known it is difficult to derive the effective Lagrangian from QCD for the general cases, and one 
popular way is to construct the lowest order operators in the Lagrangian by using the global symmetries 
of QCD as a guiding principle [l4T |. 
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It is well-known that pseudoscalars play a special role in the low energy hadronic QCD physics: they 
are the Goldstone particles which tightly relate with the spontaneous chiral symmetry breaking. In 
QCD quark model with three flavors, it is observed that the quark fields possess a U(3)ji x U(3)l chiral 
symmetry which is larger than the SU(3)r x SU(3)l- However, if the chiral symmetry breaking is 
a correct concept to the pseudoscalar sector, it is expected that the lowest Goldstone particles should 
form a nonet for the symmetry breaking of U(3)r x U(3)l — + U(3)d, instead of an octet observed in 
experiments which is consistent with the chiral symmetry breaking SU(3)r x SU(3)l — * SU(3)d- This is 
the famous Ua (1) problem for QCD quark model. A workable model should describe not only the chiral 
symmetry breaking, but also the solution of {/a(1). 

One of the early efforts is to extend the large N c effective chiral Lagrangian [ID, EB] proposed to solve 
the Ua(1) problem by including the sector for the glueball candidates, like those authors of [13, EH EH 
did by coupling the scalar glueball field to the logarithmic determinant of chiral scalar field. Another way 
is to extend the chiral Lagrangian of meson states by including a glueball field, as shown in references 
[20| . The drawbacks of these models lie in the fact that most terms in the models are put by hand or 
QCD global symmetries arc not realized in a transparent way. A model motived from QCD might be a 
better one. 

An alternative chiral symmetry model, proposed by 't Hooft [2lj j . can solve the ?7a(1) problem by using 
the instanton effects [22| . In response to the criticism from [23| , 't Hooft showed that indeed by including 
a spurion field (which will be interpreted as glueball state in this paper) coupling to the determinant of 
chiral scalar field, not the logarithmic one, the effective Lagrangian is consistent with the global chiral 
symmetry and the ?7a(1) symmetry of QCD. Meanwhile, in this form, the solution of the Ua{1) problem 
can be realised in a spontaneous way, although it is not necessary to be. The Ua(1) problem is neatly 
solved in this model. The rf corresponds to the rf observed in experiments. The determinant term 
induced from instanton effects can split the octet and singlet of the nonet and can make rf heavier than 
rf . Therefore this effective Lagrangian is well justified, and we will use this effective Lagrangian as a 
framework to describe the mixing and interactions of pseudoscalar glueballs and r\ mesons. 

In this work, we study the vacuum structure and the pseudoscalar spectra (77, rf , and psuedo-scalar 
glueball states) and their mixing. By making a connection between the FKS formula and the effective 
model, we find some correspondences between the effective parameters and the quantities computable in 
lattice QCD. We also investigate the theoretical uncertainties and the lattice computations of the model 
parameters. We explore the parameter space of the model from the constraints of the ++ and h sector 
and study some phenomenologies related to the pseudoscalar glueball candidates. It is found that three 
r\ states, i.e. 77(1295), 77(1405), 77(1475), can be glueball candidates if we assume that the lightest ++ 
glueball has a mass 1710 MeV. While there is no v glueball candidate if we assume that the lightest 
++ glueball has a mass 660 MeV. 

The paper will be organized as follows. The 4D linear chiral Lagrangian with a term induced by 
instanton effect is introduced in Section II. The mass matrices for the ++ and h sector are given in 
Section III. The theoretical uncertainties and lattice results in determining k are explored in Section IV. 
The effective mass of gluon is also explored in this Section. We analyze the parameter space determined 
from vacuum conditions and ++ and h scenarios in Section V. The phenomenological consequence of 
pseudoscalar is also studied in this Section. We end this paper with discussions and conclusions. 

II. PSEUDO-SCALAR GLUEBALL IN THE 4D LINEAR/NONLINEAR a MODEL 

The linear a model has been utilized to accommodate the chiral symmetry breaking in QCD since 
1960's [24]. In order to accommodate the glueball candidates, we must extend the ordinary model by 
including a glueball field. The effective Lagrangian is provided in [2l[ and can be formulated into two 
parts, the first part is symmetric and the other part includes explicit symmetry breaking terms. The 
effective Lagrangian reads as 



C — Cs + £-b 



(2.1) 
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The symmetry part can be put as 

C s = Tr[dX^ • 8X] - m 2 x Tr[X* X] + dY* ■ dY - m Y Y* Y (2.2) 
-*™-Tr[X*X]Y*Y - ^L T r[X^XX^X] - ^LTr[X^X] 2 



_±v_(Y* Y f-(^Ydet(X) + h.c)j+---, 

where X is a 3 x 3 bi-fundamental field which is to describe the scalar mesons (both the ++ nonet and 
the 0~ + nonet) and Y is a complex scalar (also called "spurion field") which is supposed to describe the 
pure gluon-bound states here. 

The terms which explicitly break symmetry breaking are given as 

C B = {^Tr[B ■ X] + h.c)j + Y + h.c)j + • • • , (2.3) 

which are introduced to describe the mass pattern of Goldstone particles due to the explicit chiral 
symmetry breaking. In the chiral Lagrangian, the parameter B is a matrix which is proportional to 
the quark mass and is parameterized as B = Bovxmq, where m q is a mass matrix of quarks. Similar 
to the B parameter, the D parameter should be proportional to the effective gluon mass and can be 
parameterized as D = D a ^-m e g ^ . 

It is helpful to count the mass dimension of each variable and quantity in the Lagrangian. The fields 
X and Y have mass dimension [GeV]. The parameters B and D have mass dimension [GeV] 3 . The 
parameter k is dimensionless. In this paper, for the sake of simplicity, we assume these parameters are 
real. 

The Lagrangian £$ is invariant under the chiral flavor symmetry SU(3)l x SU(3)r transformation: 

X ^ X' = U R X U f L , (2.4) 

where Ur and Ul are group elements of the SU(3)l x SU(3)r symmetry, which can be parameterized as 
exp(i9 R T a ) and exp(i9 L T a ) , respectively. T a with a = 1, • • • ,8 are generators of SU(3) Lie group, and 6 a 
are parameters of the group space. The Lagrangian Ls is also invariant under the Ua(X) transformation: 

X exp(i6 A T a )X, (2.5) 

Y -» Y 1 = exp(-i 3 A )Y . 

Where T° is the 3x3 identity matrix and 9 a is a phase angle. 

The chiral symmetry SU(3)l x SU(3)r of the Cs is broken by the nonvanishing vacuum expectation 
value vx of X field to the diagonal SU(3)d- Meanwhile, the Ua(1) symmetry is also broken by the 
nonvanishing vacuum expectation value v Y of Y field. Therefore, at low energy region, the relevant 
degree of freedom is Goldstone bosons. Therefore we can parameterize relevant Goldstone fields in a 
nonlinear way as 

X -» v x exp(i—T a ), (2.6) 

Vx 

Y -» ^-exp(i — a), 

3 Vy 

where i] a denote the Goldstone particles which break the SU(3)l x SU(3)r to SU(3)d, a denotes the 
Goldstone particle which breaks the Ua(X) symmetry. 

The vacuum expectation values in the parameterization of Goldstone field have two functions: 1) make 
the field rf and a having the right dimension in 4D; 2) make the kinetic term of all field being normalized 
to one. The advantage of the nonlinear parameterization is that the quartic terms in the Higgs potential 
do not affect the Goldstone interactions and the model can be easily compared with the chiral Lagrangian 
method. 
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We will concentrate on the study of the effects of the instanton induced determinant term, and hereafter 
we drop the term proportional to X XY , which also describes the interaction between mesons and glueballs. 
For the sake of convenience, here the vacuum is assumed to have the following form 

(X) = V X /3X3 , (2.7) 
(Y) = ^. (2.8) 

Here we already assume that v Y is real and has no complex angle. 

To have the spontaneous symmetry breaking, it is required that the vacuum should be stable. Therefore 
the parameters in the potential should satisfy the following two constraints 

k 

- 3B + 6m 2 x v x + -v Y v 2 x + 3X x v 3 x = 0, (2.9) 
o 

-D + \m\v Y + ^v 3 x + ^=X Y v Y = 0, (2.10) 



3 Y 3 A 27 



where Ax = Af + Xq 



x 



III. MASS MATRICES OF 0++ AND + SECTORS 



There are 8 free parameters in the model. In order to reduce the number of free parameters and to 
simplify the analysis, we can use the information from the ++ sector to constrain the parameter space. 

For the sake of simplicity, here we only consider the mixing between the pure ++ glueball (the modulus 
of complex Y field minus the vacuum expectation value v Y ) and the singlet <r° field. Then the mass matrix 
of ++ scalars, which determines the masses of glueball and qq states, is given as 



M h = I I ■ (3-1) 



M\y Myy 



While the mass matrix elements are given as 



M xx = 6m 2 x + 9X x v x + hv Y v x , (3.2) 

M YY = 2m Y + ^-v Y , (3.3) 

M XY = kv\. (3.4) 

We can use the physics mass of ++ glueball and qq states and the mixing angle to represent these 
elements as 

M 2 q cos 2 a + M 2 g sin 2 a = M xx , (3.5) 

M 2 q sin 2 a + M 2 g cos 2 a = M YY , (3.6) 

\ (K - K) sm 2a = M xy ■ (3-7) 



Then by using Eqs. (12. 9112.101 and Eqs. (I3.5H3.7|) with six inputs ( the experimental determined values of 
Vx, the mass of ++ glueball candidate M gg , the mass of qq meson state M q g, the assumed mixing angle 
a, the assumed value of v Y and Tn?J* ), we can solve out the variables m x , m Y , X Xl Ay, and k. 

In order to analyze the symmetry breaking pattern, we choose two scenarios to scan the parameter 
space of Ax and Ay . The masses of scalars are chosen as given in Table |TTJ Scenario one is motivated 
from the possibility that <r(660) might be a glueball rich state (2f| [2(|. Scenario two is motivated from 
that /o(1710) is a good glueball candidate, as proposed by 27] and supported by the lattice computation 



G 





0++ 


0++ 


Scenario 1 


0.66 


0.98 


Scenario 2 


1.71 


1.37 



TABLE II: Two scenarios in ++ sector are considered for the study of symmetry breaking pattern and for the 
constraint on the parameter space. 



[40] . It is found this scenario is also favored in the recent study under the context of AdS/QCD [28j| , We 
assume the mixing angle is a free parameter for these two cases. By requiring that the potential should 
have a bottom when X — > oo and Y — > oo, we can impose the following two conditions Xx > and 
Ay > 0. At the same time, we require that the value of Ax and Ay should not be too large to violate the 
unitary conditions then we impose two conditions Xx < 40 and Ay < 40. 

From the solutions to the five equations, we observe that there are four patterns of symmetry breaking in 
the parameter space: 1) Symmetries are broken by the conditions m x < and m Y < with k > 0, which 
corresponds to the case where these symmetries are broken independently to each other. 2) Symmetries 
are broken by the conditions m x < and m\ > with k > 0, which corresponds to the case where 
the chiral symmetry breaking triggers the J7a(1) symmetry breaking. 3) Symmetries are broken by the 
conditions m x > and m Y < with k > 0, which corresponds to the case where the Ua(1) symmetry 
breaking triggers the chiral symmetry breaking. 4) Symmetries are broken by the conditions m x > 
and m Y > with k < 0, which corresponds to the case where instanton effects trigger both the chiral 
symmetry breaking and the J/a(1) symmetry breaking. The fourth one is a new pattern which can be 
related to the study on the chiral symmetry breaking by instantons [2^ |. 

To analyze the ^ sector, we use the nonlinear realization formula. By substituting the nonlinear 
parameterization of Eq. (|2.T[) into the Lagragian of Eq. (|2.1[) . we obtain the kinematic terms of 7r°, r]°, 
rj 8 , and a fields, which determine their masses and how these pseudoscalar mix with eath other. Up to 
quadratic terms, the Lagrangian reads as 

C = -dir° ■ dir° + ^dri q ■ drj q + ^dri s ■ drj s + ^da ■ da (3.8) 

i 0\2 -"9 I \2 B s / N 2 ^ D . . 2 

- ttH 71 " ) ~7T-(Vq) -ir-iVs) -7, — («) " 

ZV q ZV q ZV S IVY 

where we replace 77° and 77 s by the combination of r\ q and rj s 




vx V 3 v q V 3 v s 

— - Vs a -VI-- (3 - 10) 

vx V 3 v q V 3 v s 

Obviously, neither tt° nor rf mix with a and rf . Only the last term, originating from the instanton 
effects, leads to a mix between rf and the ^ glucball a. Such a mixing term splits the mass of rf from 
that of the rest of pseudo-Goldstone bosons and provides a resolution to the Ua{^) problem of QCD. 
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From the Lagrangian, it is straightforward to read out the mass matrix as 



M = 
gqs 



gg 



kv 2 ^ 
g vy 

^kVqVs 



3 kv g 



' 2 

I 2 + 2 - 

qq ' 3 

V2 



3 kv q v s 



3 kv q 



\ 



kvyV s ^kvyVg 



kv Y v Q m? + Uv 2 ^- 



(3.11) 



The m 2 is equal to 3D/vy. In order to reduce the number of free parameters, we will assume that 



Bn 



Dq- We take the mi = B/v q = Bo(m u + nid) = m 2 , and mj s 



2M 2 



Ml 



It is instructive to compare the mass matrix determined by the effective field theory method with that 
in the FKS formalism [30(. Here we only consider the leading order appromixation. The terms in the 
first 2x2 components of the mass matrix can be compared directly to the FKS formalism ( We provide 
some crucial steps to derive the correspondence in Appendix ). We can have the following one-to-one 
correspondence: 



qq 

„2 



-kvyV, 



-kv 2 ^- 
3^ v s 









0|m u ui75 u + 


v q 








\m s sij 5 s\ri s }; 


v s 












-iff 

V s 


\^GG\ % ). 



mdd«7 5 d|?7 q ); 



l^GGM; 



It is straightforward to have the following relations when FKS formalism is extended to include 
doscalar glue ball 



in,, 



kv 2 



Vy 



kv q V s 



V5 

Vy 

V2 



4,n 



GG\a) 



(0\d^ q \a 



Vy 47T 



1 



V 8 ' ' ' ""' ' Vy ' 47T 

Then we can have the following solutions in the leading approximation 



ii>7, 



, a. 



(0\^-GG\ Vq ) 
{0\^GG\ Vs ) 

(0\^GG\a) 



0, 
V2 



kvyVqV 



q"s , 



■^kVYV q , 

(0\d^ q \a) = (0|^J 5 » 



kV 2 Ve 

3 9 



3.12) 
3.13) 
3.14) 

3.15) 
3.16) 

pseu- 

3.17) 
3.18) 
3.19) 

3.20) 
3.21) 
3.22) 

3.23) 



Thanks to the relation given in Eq. (|3 . 23[1 . it is possible to extract the value of k from the lattice 
calculation. Due to the relation between the effective field theory and the FKS formula, it is also possible 
to extract the value of k from the pseudoscalar meson spectra. In the following section, we will discuss 
the computation of k in QCD by considering the instanton effects and its uncertainties. 
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IV. THE k PARAMETER AND THE EFFECTIVE GLUON MASS m% ss 

The generating functional induced by one instanton can be deduced from QCD, as shown by 't Hooft 
[22| . which reads 

N f / - \ 
k t e e Aet{i>[{x)^ R {x)) , (4.1) 

where ip denotes the quark fields, k t is a constant, x is the position of instanton, the / denotes the index 
of flavor space, and 9 is the vacuum angle. Following the de-clustering assumption and summing over all 
the contributions of both instantons and anti-instantons, we can arrive at the effective interaction which 
reads as 

exp | J k t e e det (^p f L (x)ip f R (x)j + /i.e. j . (4.2) 
To relate this effective interaction with that of the meson fields, we need redefine the following variable 



A 2 

match 



(4.3) 



where A matc h is the matching scale between the fundamental theory and the effective theory, which is 
introduced from the dimensional analysis. Then we can express the effective interaction of mesons induced 
from instantons as 

exp | jf k t AZ{ ch e e det (cf>(x)) + ft.cj . (4.4) 

In order to formulate the solution of the Ua(X) problem within the framework of spontaneous symmetry 
breaking, it is instructive to introduce a field Y to realize the continuous symmetry plj . Therefore the 
effective interaction arrives at the form 

exp | J ^Ydet (<f>(x)) + fe.c.j , (4.5) 
which is given in Eq. (|2.2[) . where the effective coupling k is expressed by k t as 



fi£vA 2JV/ 

k = 1 match . (4.6) 



vy 

Implicitly, this Y field is related with the instanton configurations of QCD. Compared with other forms 
of effective interaction proposed in literatures, the advantage of this interaction is that the anomalous 
Ward identities can be easily satisfied. 

Now we address the strength of the effective coupling k t . k t can be obtained from the one-loop 
calculation. For the SU(Nc) group, the contribution of gluons can be formulated as [3lJ 



2N C 

7T 2 (iV c -l)!(AT c -2)! \g(iif 

where N c is the number of color. Obviously, in the large N c limit, k t vanishes. The contribution of Nf 
chiral quarks can be formulated as [13] 

k q = 2 3N fn 2N f p 3N fe -%N f ln(pp)-2N fa (±) _ 
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Combining the contributions of gluons and quarks (the integration of the product of kf and k\ over 
instanton size p), we arrive at 



h = Ci 
ft = 4 



n 3N l~ 5 



4tt 2 



2 \ 2^Vc 



5(M) 2 , 

e -a(l)-2(JV e -2)a(i)+2iV / Q((i) 

(iV c - l)!(iV c - 2)! 



To obtain this expression, we have utilized the renormalization group equation of gauge coupling 



Sir 2 



8tt 2 



<?(i/p) 2 g{n? 



li 



ln(np) . 



(4.9) 
(4.10) 

(4.11) 



While the constants a(l) and a(h) are dependent on the regularization scheme and are given in pH ]. 
According to the study of Geshkenbein and Ioffe in 32], it is better to use the results of dimensional 
regularization. Then these two constants are given as 



1 16 

ail) = 8R+-ln2 

V ' 3 9 

,1, 1 17 

q(-) = 2R - -ln2 - — 



2' 6 72 ' 



12 

0.248754 . 



2tt 2 ^ 

2 



(4.12) 
(4.13) 

(4.14) 



Below we address the uncertainty of kt . There are a couple of theoretical uncertainties in determining 
the value of k t : 1) the first uncertainty is from the regularization scheme. For instance, the value of 
k t in dimensional regularization scheme can be one thousand times larger than its value in Pauli- Vilas 
regularization scheme. 2) The second uncertainty is from the infrared cutoff in the size of instantons. 

In Table we tabulate the dependence of kt on the infrared cutoff pm of the size of instantons. 

In order to explore the dependence of kt on the energy scale, we divide the integration of Eq. (|4.9p into 
three interval integration: 

i 

k\ = [ mb f K (N c ,N f )dp, (4.15) 
J o 

k\ = J™* f K (N c ,N f )dp, (4.16) 

fc ( 3 = J''" f K {N c ,N f )dp, (4.17) 

where mf, and m c are the masses of b and c quarks. The final value of k t is the sum of these three parts. 
The integrand function Jk{N c , Nf) can be read out from Eq. (|4.9p . In these three integration intervals, 
the effective active degrees of quark flavors, Nf, are different. 

There are several comments on the Table pil[) in order. 1) One fact is that the dominant contribution 
to kt is from the infrared regions as demonstrated by the value of k\ and k\. 2) The value of k t is very 
sensitive to the change of pm, for instance, when pm changes from 1 to 2, the value of kf increases by 
10 3 times, as shown in the 4th and 5th columns. 3) Values of kt in the dimensional regularization scheme 
are larger than those in the Pauli- Villas regularization scheme by a factor e 8 ' 8 . 

The value of k further suffers an uncertainty from the determination of AJ^^ and vy- Below we will 
assume that A matc h = 1 GeV. 

Table linl explicitly shows that it is difficult to determine the value of k t and k from pure theoretical side 
due to our ignorance on the value of pm. If we believe that the dimensinal regularization is the correct 
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Contri to kt 






PlR = 1 


Pir = 2 


k t 

PIR = 1 


k t 

Pir = 2 


Dim sheme 


1.30 x KT 4 


1.07 x 10 3 


7.45 x 10 1 


1.11 x 10 6 


1.15 x 10 3 


1.11 x 10 6 


P-V scheme 


1.97 x 10" 8 


1.61 x 10" 1 


1.12 x 10" 2 


1.68 x 10 2 


1.73 x 10" 1 


1.68 x 10 2 



TABLE III: The values of kt in different pirs and regularization schemes are shown. The Dim and P-V schemes 
denote for the dimensional regularization scheme and the Pauli- Vilas regularization scheme, respectively. The 
Pir has the unit as GeV . 



scheme and to have k ~ 60, it seems that the pm should be around 0.5 GeV^ 1 or so, which corresponds 
to Am w 2 GeV. Here Am is the matching scale of QCD and the effective chiral Lagrangian given in Eq. 
(|2.ip . While in the Pauli- Villas regularization, the predicted k is too small to solve the J/a(1) problem. 

From now on, we address the effective gluon mass m e g ^ . It is well known that in the perturbation and 
large momentum region, the mass of gluon vanishes. In the infrared region and large distance region, the 
effective mass of gluon might be non-vanishing. Currently there are two methods to extract the effective 
mass of gluon: one is from the Dyson-Schwinger equations as demonstrated in references [33l l34l ]. the 
other is from the lattice calculation (as demonstrated in references [35l l36l. l37j). 

From the method of the Dyson-Schwinger equations, as shown in references [H, 34], the gluon propa- 
gator can be computed and the effective mass of gluon can be extracted. Depending on the ansatz of the 
three-gluon vertex, two types of solutions for the effective mass of gluon can be derived. The first type 
of solutions can be parameterized as 



m 



(« 2 ) = WM^^VM^)]- 1 - , (4.18) 



. g j J »m L ~v A2 ,1 -V A2 

Lass {my* 

The second type of solutions reads as 



then the effective squared mass (m^-^) 2 (q 2 ) has a logarithmic running behavior, while 71 = |(1+C2 — c\). 



I em 2 /2N "V n f q 2 +pm 2 f pm 2 ., ^ in , 

K JJ ) (q ) = q2 + mo2 M — j2 — VM^HP . ( 4 - 19 ) 

then the effective squared mass (m^^) 2 (g 2 ) has a power law running behavior, while 72 = | + ^f 1 - The 
A is the mass scale of QCD, which is taken as A = 0.30 GeV. 

We show these two types of solutions in Fig. ©. In Fig. ([2k), the parameters are taken as m\ = 0.3 
GeV 2 , p = 1.007, c\ = 0.15, and C2 = —0.9635. In Fig. (J2J)), the parameters are taken as m§ = 0.5 
GeV 2 , p = 1.046, ci = 1.1, and c 2 = -1.121. 

From Fig. ([2]) and the analysis given in reference [34 |. we can read out that the effective mass of gluon 
can expand in a range 0.5 GeV < ra e J^ < 1.2 GeV when q 2 — > 0. 

The lattice method can also extract the effective mass of gluon. For instance, the references [3?], [H[ 
present results where the gluon mass in SU(3) case is 0.3 GeV ~ 0.5 GeV or so. The results are 
claimed to be consistent with the so-called decoupling solutions found for Dyson-Schwinger and functional 
renormalization group equations. Meanwhile, as shown in 39], it was found that there is very little N c - 
dependence in the gluon's propagator over the whole momentum range, so is the effective mass of gluon. 

Combining the results given by these two computational methods, in this paper, we will assume that 
the effective mass of gluon is in the range 0.3 GeV < m e J^ < 1.2 GeV. 



11 




f[Gef; 



(a) 



(b) 



FIG. 2: 27ie dependence of the effective gluon mass m 2 on the external momentum q 2 is shown. Plot (a) depicts 
the case where m 2 has the power law running behavior, and plot (b) depicts the case where m 2 has the logarithmic 
running behavior. 



V. NUMERICAL ANALYSIS 



To do numerical analysis, we will take the following inputs: 

v q = 0.092GeV , v s = V2v q , m„ = 0.003GeV , m d = 0.007GeV , 
m s = 0.15GeV , m„ = 0.1395GeV , m K = 0.493GeV 



(5.1) 



Below we concentrate on the analysis of the allowed parameter space of my* , k, and k vy. First 
we study how the ++ can constrain the parameter space. To scan the parameter space, we allow the 
effective gluon mass to vary in the range 0.3 GeV < J7i|" < 1.2 GeV, the vacuum expectation value vy 
to vary in the range 0.001 GeV < vy < 1 GeV for the first scenario and 0.001 GeV < vy < 1.5 GeV for 
the second scenario, and the mixing angle to vary in the range — ? < a < ^. 

Fig. (0 shows the parameter space allowed by the vacuum constraints, untarity bounds, and the ++ 
inputs. Comparing [3(a) and[3jb), we observe that the allowed region for nfj* and kvy is larger in the 
second scenario than that in the first scenario. It is similar for k, as demonstrated by [3(c) and [3(d). One 
interesting observation is that the value of k can be either positive or negative, which can realize the 
fourth pattern of the chiral symmetry and £7a(1) symmetry breaking. 

Now we study the parameter space allowed by the ^ sector and treat k, vy, and fny' as parameters 
with no bounds. Substituting the inputs given in Eq. (|5.ip into the mass matrix of pseudoscalar, we can 
arrive at the following mass matrix which describes the mixing of pure v glueball state, qq state, and 
ss state: 



M 



1.0 x 10" 3 -^ -+ 
9.8 x 10 



1.95 m e / f 
3 k 

4.89 x 10~ 3 k 



9.8 x 10~ 3 fc 
1.95 x 10~ 2 + 8.67x 1Q- 2 kv Y 
4.34 x 10~ 2 kvy 



0.47 



4.9 x 10~ 3 fc 
4.34 x 10^2 kvy 
2.17 x 10~ 2 kvy 



.(5.2) 



In order to make the solution of the ?7a(1) problem explicitly, it is also useful to rewrite the mass 
matrix in the basis of a, 770, and 77s- The matrix reads as 



1.0 x 



M 2 = 



10-3 Jl 

vy 



1.0 X 

1.6 x 



I- 1.95 m e /f 

1-2; 



io- 2 fc 

10" 3 fc 



1.0 x 10" 2 fc 
0.168 + 0.1 kvy 
-0.21 + 1.6 x IQ- 2 kvy 0.32 + 2.5 



1.6 x 10~ 3 fc 
-0.21 + 1.6 x 10~ 2 kvy 

X IQ- 3 kvy 



(5.3) 



This numerical mass matrix reveals several transparant but interesting facts on the solution of the J7a(1) 
problem, which will be listed in order below. 1) Comparing the 22 and 33 element, we can know that 
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£■ v =(0.75, 1 .DO 1 Ge 

* v ,—(0.50, 0.75 ] Ge\ 
v,= ;0.25, 0.50 lGe\ 

• v„=;0.00, 0.25 JGe\ 



C= v,= (1.U, 1.50] GeV 
* v,= (0.75. 1.13 1 GeV 
" ■ v,= (0.38, 0.75 ] GeV 
v,= (O.0b, 0.58 ] GeV 



. ■ -n "... <™ * ■' ■ 



C 0.25 0.5 



nfffGeV) 



(b) 





</,= (□. 75, 1,00 1 GeV 
O v,= (0.50. 0.75 1 GeV 
fl= v,= (0.25, 0,50 1 GeV 
* v,= (0,05. 0.25 ] GeV 



nff(GeV) 



(c) 



:\;<^.;.V>;.--' ' 



v,= (0.375, 0.75 ] Ge\ 



3 , 



,-(0.05, 0.38 ] GeV 




nff(GeV) 



(d) 



FIG. 3: The parameter regions from the constraints of vacuum conditions and the ++ sector are shown. Plots 
(a) and (c) are devoted for the scenario one, and plots (b) and (d) for the scenario two. 



if k = 0, then m^o is always smaller than m„s. 2) Generally speaking, if k << 10 3 , the mixing of pure 
glueball a and rf is negligibly small, which indeed supports the two angle scheme to describe the mixing 
among pure glueball state a, r\ q q and riss, as demonstrated by 12 and 13 elements. 3) In order to solve the 
C/a(1) problem, the product of k and vy must be large enough (say 8 ~ 13) to lift the m^o and realize the 
scenario with m v o > m v s, as indicated by the 22 element. While due to the suppression factor 2.5 x 10~ 3 
in the 33 clement, the mass of 77 s is not sensitive to the change of k vy. 4) From the 23 element, we can 
read that there is a minimum mixing scenario for rf and 77 s , which correspond to kvy 13. 5) If 
is very large, the glueball state might decouple from the solution of ?7a(1) problem. 

In order to examine how the parameter k controls the spectra and mixing, we select the following two 
bench mark cases to study the dependence of mass spectra and mixing on parameters. These two bench 
mark cases assume that 77(1405) and 7/(2190) are glueball candidates, respectively. The values of vy and 
THgff are selected from the solutions by using the least x 2 method. Such solutions are possible since 
the mass matrix of pseudoscalar includes three free parameters while we can require that the eigenvalues 
of the mass matrix correspond to the masses of 77(548), t/(958), and the mass of the physics glueball 
candidate, which is equivalent to impose three conditions. Table IIVI lists the values of these two bench 
make cases. 

In Fig. we use a pie with three colors to denote the three components of each spectrum: red, 
blue, and green denote glue ball, qq (or 77 ), and ss (or 77 s ), respectively. The fraction of these three 
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Vy 


m% ss 


k 


Bench case 1 7/(1405) 


0.234 


0.824 


37.75 


Bench case 2 7/(2190) 


0.133 


2.160 


58.72 



TABLE IV: Two bench mark cases to show how the spectra of pseudoscalars and their fraction are affected by 
the parameters. The best fit parameters are determined by using the least ~xf method. We select only one solution 
for each bench case, though the solution for each case is not unique. 



Mass Spectra (GeV) 
3.0 h 



Mass Spectra (GeV) 
3.0 - 



ta(k) 



(a) 



Mass Spectra (GeV) 
3.0 h 
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(b) 

Mass Spectra (GeV) 



#••••••••#• 



••••• 



•••tit 



ta(k) 



(c) 



(d) 



FIG. 4: TTie dependence of the spectra and their fractions on parameter k is shown. The y axis is the mass of the 
spactra, and the x axis is the parameter k. Each pie in the plot denotes the component of the glueball, qq ( or rf 
), and ss ( or rf ) states. The figures (a) and (b) are in the basis of a, qq, and ss. The figure (c) and (d) are 
in the basis of a, rf , and rf . The red part of each pie denotes the glueball component. For figures (a) and (b), 
the green part denotes the qq component, and the blue one denotes the ss component. For figures (c) and (d), 
the green part denotes the rf component, and the blue one denotes the rf component. The y axis is deliberately 



chosen to located in the best fitted k for each plot. 
of 7/(548), 7/(958), and 7/(1405) (or 7/(2190);. 



The red, blue and green lines represent the experimental values 



components is determined by the eigenvector of each eigenstate. For instance, if the eigenstate of ?/(548) 
is found to be |?j(548)) = V ng \gg} + V m \qq) + Kj S |ss), where V V i,i — g,q,s is defined as the component 
of eigenvector. These components must satisfy the normalization condition as V^, t = 1. The fraction 
of each component in the pie is determined by V^ i: Fig. ([4^) and Fig. §4)p) depict the bench case for 
rrir, = 1.405 GeV, and Fig. (do) and Fig. gH) for m v = 2.190 GeV. 
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From these plots, we can read out the fact that when k is small (say k < 15 or so), the mixing is 
negligibly small and the mass eigenstates are almost pure states. Then the state with dominant qq (or 
?y ) component is lighter than that with dominant ss (or 77 s ). When k increases (say 15 < k < 60) , the 
mass of qq (or 77°) state increases quickly, while the mass of ss (or 77 s ) state is not sensitive to the change 
of k. When k is large enough, the mixing between qq and ss becomes significant and can realize the 
masses of 77(548) (which is dominant rf component) and rj 1 (958). If k increases further (say k > 60), the 
lightest state (77 s ) is almost insensitive to the change of k, while the masses of the heaviest and second 
heaviest states become larger. It is also noticed that when k > 60, the heaviest state can be either 
glueball dominant state or the qq dominant state. 

We also observe that in the solutions, the 77' (958) contains more glueball component in the bench case 
one than in the bench case two. The bench case two demonstrates the decoupling limit of glueball state 
in the solution of the Ua(1) problem. This fact indicates that if glueball candidate is light (say less than 
1500 MeV), ?y'(958) might has considerable amount of glueball component. While if glueball candidate 
is heavy (say heavy than 2000 MeV as predicted from lattice computation), then 77' (958) might has no 
significant amount of glueball component. However, such an amount of glueball component in 77' (958) 
depends on the parameters of the solutions. 

For the solutions, we observe one interesting fact: in order to solve the J7a(1) problem, the heaviest 
state is always the glueball dominant state, and the second heaviest state is always rf dominant state. 
There is no solution for the case where 77' is dominantly glueball component. 

From these two bench cases, we know how the C^i(l) problem is solved and the role played by the 
glueball candidates. Here and below, we extend to study the region of parameter for nine possible 
psedudo-scalar glueball candidates, as given in Table [U 

Fig. © is devoted to show the solutions in the parameter space for each case. Fig. ([5]) demonstrates 
several facts: 1) Although the parameter k can extend from 20 to 180 for these nine cases, the product 
of k and vy is situated in a quite narrow range (from 7 to 18), as shown in plot (a). Generally speaking, 
the small vy corresponds to the large k, vice versa. This indicates that the term k vy in the mass matrix 
given in Eq. (|5.3[) has more influence than the term in determining the solutions. By using the relation 

between the effective field theory and the FKS method, kvy is related with the values of (0| GG\r] q ) and 
(0| GG I r/s). Then from Fig. (|5]c), we can read out these transition elements: (0| GG\rj s ) is within 
(0.02,0.06) GeV 3 and (0|££ GG\r] q ) is within (0.04,0.12) GeV 3 , which is simply twice of (0|fj GG\r) s ). 
These results are reasonable when compared with lattice results [13,|4l|. 2) When the product of k and 
vy is fixed, heavier glueball candidates have solutions with larger value of rn e Jf. Roughly speaking, the 
correlation between the product of k and vy and m e Jf is linear, quite similar to that between k and 
m e Jf. The value of k is directly related with the transition element (0|^ GG\a), which is determined to 
be within the range (0.01,0.1) GeV 3 , as demonstrated by ((SJl). 3) A heavier glueball has a wider range 
of k and k vy . 

Solutions in Fig. O demonstrates a lower bound for m e g ^ with vn?J* > 0.5 GeV, which is larger than 
the lower bound determined from lattice and Dyson-Schwinger methods. It is also worthy of mentioning 
that the information on the effective gluon mass can be helpful to rule out lots of solutions. For instance, 
solutions with Trig** > 1.2 should be excluded if we treat the upper bound of Tn e J' from lattice and 
Dyson-Schwinger methods seriously. 

Below we calculate the branching fraction of pseudoscalar glueball candidates (which is denoted as G 
here) decaying to two photons, decaying to two electrons, and to two muons. We assume that glueball 
only decays into these channel via its quark component [42]. With this assumption, the ratio of the 
G — > 77 width over the 7r° — > 77 one is expressed as 



For the G — > £ + £ decays, it is possible to determine their widths by using the available ir° — > e + e and 




(5.4) 
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(c) (d) 

FIG. 5: The solutions for nine possible pseudoscalar glueball are demonstrated. The x-axis is for the effective glue 
mass. The y-axis is the product of k and vy for (a), k for (b), (0||j GG\rj s ) for (c), and (0| GG\a) for (d). We 
only show solutions with vy < 0.4 GeV. The convention for marker in (c) and (d) is the same as in (a) and (b). 
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(5.5) 

glueball candidates, m gg and the corresponding 



These decay widths only depend on the mass of 
eigenvectors. 

Fig. © is devoted to study these decay widths, which might be useful to identify glueball candidates 
from future experiments. We notice that T(rj — > 77) can change in the range (0.36, 0.45) KeV (slightly 
lower than its experimental value 0.51 KeV) [l[ and Tin' — > 77) can change in the range (3, 5.5) KeV( 
its experimental value 4.28 KeV is within this range [1J ), both of which are not very sensitive to the 
change of k. However, as we can still read out from Fig. ([6^,), for each case, when k is large, T(rj' — > 77) 
decreases slightly. The underlying reason is that when k increase, more glueball component enters into 
77', which reduces the quark components in effect. 

While in each case, T{G — > 77), T(G — ► e + e~), and T(G — > /i + /i _ ) increase with the increase of k, as 



1G 





(c) (d) 

FIG. 6: The dependence of decay widths of glueball candidates on parameter k are shown. In order to compare 
and contrast, we also show the decay width of rj' to two photons. In plot (a), the convention of markers is the 
same as the rest of plots. 



demonstrated in Fig. (JSh-Oi). These results demonstrate that it is quite a challenge to measure these 
processes from experiments. For each case, the relation between decay width and the parameter k is 
simply linear, since large k means that more quark states mix into the glueball state, as revealed by the 
mass matrixes Eqs. (I3.11I5.2I5.3|) . Another remarkable fact is that the heavier the glueball candidate, 
the wider its decay widths might expand in the model. 

It is also interesting to study how the mass splitting between the ++ and h glueball. Fig. ([7]) 
is devoted to study the mass splitting between the ++ and h glueball states. We define Am g — 
m ++ — m -+. In the ++ sector, the mass of glueball candidate is fixed by assumption. In the h 
sector, we identify the eigenstate with the largest V?,i — i],r)',G as the glueball candidate. From Fig. 
([7]), we observe that in both cases, the mass splitting is simply controlled by the effective gluon mass. 
In the first case, the larger effective mass of gluon corresponds to a larger mas s sp litting and the h 
glueball state is heavier than the ++ state which is similar to the lattice results [4(| , while in the second 
case, the situation is opposite. Therefore, if the mass of the glueball candidate is large than 1295 MeV, 
it is difficult to realize in the first scenario. While for the second scenario, only those states with mass 
smaller than 1760 MeV can be realized. Another interesting fact as shown in Fig. ([7]) is that the mass 
splitting is almost independent on the parameter k. 
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(a) (b) 
FIG. 7: Mass splitting between the scalar and pseudoscalar glueballs is demonstrated for two ++ scenarios. 



VI. DISCUSSIONS AND CONCLUSIONS 



In this work, we use the effective Lagrangian derived from the instanton effects to study the pseudoscalar 
glucball candidates. At the leading approximation, we observe that the results obtained in this model 
are consistent with other methods, like FKS method and the lattice method, which indicates that this 
effective Lagrangian captures the most important features of QCD. 

We notice that the ++ and h sectors favor different regions of the parameter space, therefore the 
physics regions of parameters for both sectors are constrained. Without considering the Regge excitations, 
we find that the parameter space for the first scenario given in Table Ql] for the ++ sector can not be 
consistent with the known masses of pseudoscalar glueball candidates while the second scenario for the 
++ sector can. In the second scenario given in Table [IT] for the ++ sector, the pseudoscalar glueball 
candidates cannot be heavier than 7^(1760), therefore only three candidates, i.e. r/(1295) , r/(1405) , and 
77(1475), are possible pseudoscalar glueball candidates. 

To study the parameter space of m e g ^ , we have assumed that Do = B . Currently, we do not know 
the exact relation between D and m e Jf from the underlying QCD theory. If Do is much smaller than 
£?o, then the parameter can not be constrained too much. While if assuming Do = 3 Do, we notice 

that the unitarity conditions rule out all points in the first scenario given in Table |TT] of ++ sector and 
only few points survive for the second scenario given in Table ILTl 

We have combined both ++ and h sectors to study the constraints on the parameter space. In 
the ++ sector, we only consider the mixing between the singlet <r° and the ++ glueball. In a more 
realistic treatment, we should take into account the realization of the experimental ++ spectra. Then 
it is necessary to extend the Lagrangian by including the tetraquark state [43|, [44|. Furthermore, we 
have utilized the nonlinear realization formalism to treat the ^ sector. In order to take into account 
the interaction between ++ and h sectors, maybe it's better to use the linear realization form to 
parameterize both sectors, as demonstrated by [25, 45]. However, the main conclusions given here might 
not be changed. 

We have found the correspondence between the FKS method and the effective field theory method at 
the leading order. In the effective field theory method, the mass matrix is automatically symmetric by 
requiring the Lagrangian is Hermitian. While in the FKS method beyond the leading order, the mass 
matrix is not necessarily symmetric. In order to have a symmetric mass matrix, we must invoke an 
asymmetric decay constant matrix (|8.22p . which is difficult to find its counterpart in the effective field 
theory method. 

The 4D effective field theory can be extended to accommodate higher excitations of by putting the 
masses of pseudoscalar on the Regge lines. Interactions can be introduced by invoking the assumption that 
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there exists a global symmetry for each radial excitations multiplet, as suggested in [46|. An alternative 
way is to use the AdS/QCD in the soft- wall model [47], Hfl, which can accommodate the Regge trajectories 
in a more convenient way. We will extend our analysis on pseudoscalar glueball candidates to this model 
in our future work. 
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VIII. APPENDIX 

A. The symmetric FKS formalism with glueball 

We follow the convention of Q and list the FKS formalism by including a glueball candidate. We 
extend the asymmetric mass formalism given in Q to a symmetric 3x3 mass matrix. 

The basic difference in our formalism is that we can define a vector current of gluon fields 

= ^e^A a ^d a A a p + p abc A b a A c ^ . (8.1) 

At the same time, for the quark sector, we can define two axial-vector currents, which read 

J£ 5 = 97^759, (8.2) 
Js, 5 = ^75S. (8.3) 

Therefore, beside the standard FKS decay constants (vacuum expectation values) v q and v s given below 

(QIJ&MP)) - -iv s P», (8.4) 

we have an extra definition 

(0\K"\a(P)) = - ^HyP*, (8.5) 

where P is the momentum defined in 4D. To include the OZI rule violation effects, we introduce the 
following transition elements: 



(0\K»\n q (P)) = 


fg da» 

V3 Jq 


(0\K»\ Vs (P)) = 


— 5=/» p" 


(oKskOP)) = 


— 


(o|J£ 5 KP)) = 


fi 

V2 Ja 


(o|J s M 5 K(P)> - 




WM p )) = 





(8.6) 
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The fj denotes the decay width of state j to vacuum via operator i. At leading approximation, /j 
vanishes. 

The physical states, the psuedo-scalar glueball state G, rf, and r\ are the mixtures of three basis 
ingredients a, rj q , and rj s , which can be represented as 

/ la) \ 

(8.7) 

where U is an orthogonal normalized matrix. Then relation between the decay width of physical states 
and of basis state is given as 

F = U \ tsfi f* Vq ^ 

J_ f 9 J_ f 1 




V. 



where F is defined as 



fr fr fr 
g fg 



F= % ft Y # . (8.9) 

fg fq fs 

J 7} J Tj J Tj 

Then we can have one identity and the two anomalous Ward identity, which read 

d M Ki* = ^G^G^ , (8.10) 

ry — 

d,J^ = 2im q q l5 q + -± G„„ G^ , (8.11) 
d^ 5 = 2im s s l5 s + £i G^ G» v . (8.12) 

47T 

By using (0| and physical states to sandwich these three equations, we arrive the following equation: 

(0\d^\G) (OlfyJ&IG) (Q\dMG) 

(0\d,K»W) (Q\d,r q . 5 W) (0|^^ 5 |r/> \=Ml-F, (8.13) 
(0|a M Jf*|T7) (0\dM v ) (0\dMr,) 



where M|, is given as 



Ml 



M 2 D = \ M*, ] . (8.14) 

Ml 

To obtain this equation, we just use the definition given in Eqs. (|8. 4118. 61) . 
According to Eq. (|8.8|) , we can also express right hand side of Eq. (|8. 13[) as 

/ (0\d^\G) (0\dMG) (0\dMG) \ ( (0\d^\a) (0\dMa) (0\d^ 5 \a) 

(0\d,K»W) (0\d,jl 5 W) (0\d,J^ 5 W) \=U \ (0\d,K»\r h ) (0|^A| % ) (0\dM Vq ) 
\ {Q\d M K»\r,) (0\d,J^ 5 \ V ) (0\d^ 5 \ V ) J \ {Q\d M K»\r, a ) (0|S M jJ B |»/.) (0\d,J^\ Vs ) 

(8.15) 

By using Eqs. I|8. 10H8.12"!) , the last term in the right hand side of Eq. (|8.15[) can be represented as 
{0\d M K»\a) f 3 mi qa + (0\d^\a) ^m 2 sa + {0\d^\a} \ 

{0\d M K"\ri q ) ^m^ + <0|S^| % ) + <0|9^| % ) - M 2 gqs -F', (8.16) 

(0\d^\r, s ) v s m 2 qs + (0\d^\r, s ) v s m 2 ss + (O^K^) / 
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with 



M = 
gqs 



F' = 



( £(0\d^\a) ml a + £(0\d,K»\a) 
V MO^K^s) m 2 qs + M0\d^\r, 8 ) m 2 ss + 



,2 



75 VY 



g(0\dpK»\a) 



V2 U 1 



.17) 



.18) 



The quantities mr^ are defined as 



in 



qa,qq,qs 
sa.sq,ss 



V2 



(Q\im u wy 5 u + im d d^ 5 d\a, rj q , 7y s ) , 



\tm s sj 5 s\a, r) q ,T)s) 



(8.19) 
(8.20) 



Combining Eqs. (|8.81 18.131 18. 151 18.16|) , we arrive at the symmetric mass matrix which reads as 

r 2 t-1 



where J is defined as 



/ i 11 il\ 



fg rs 
J g 1 ±q_ 



VY Vg 



1 J 



.21) 



.22) 



At the leading order approximation, J" 1 is an identity matrix. Therefore, we can match the left side 
with our mass matrix determined in the effective field theory method given in Eq. (|3 . 1 1[) . 

Obviously, in the general case, the mass matrix Mg qs is asymmetric. In order to have a symmetry mass 
matrix, the asymmetric matrix J is needed, which corresponds to take into account the higher order 
corrections. 
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